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1. Introduction 

Recently, two of the authors have considered a complex matrix model which describes 
the ensemble of branched coverings of a two-dimensional manifold [TJ. The model has been 
interpreted as a string theory invariant with respect to area-preserving diffeomorphisms of 
the target space. In this paper, we continue the investigation of this model. 

The geometrical problem we solve consists in the enumeration of the branched cover- 
ings of a two-dimensional manifold with given number of punctures. The target manifold 
is characterized by its topology, the number of punctures, and its total area. All struc- 
tures we are considering are invariant under area-preserving diffeomorphisms of the target 
manifold. We assume that the covering surfaces can have branch points located at the 
punctures. Two covering surfaces related by an area-preserving diffeomorphisn are consid- 
ered identical. 

The problem will be reformulated in terms of a lattice gauge theory of N x N complex 
matrices defined on a lattice representing a cell decomposition of the target manifold. The 
vertices of the lattice are the punctures of the target surface. The -k perturbative expansion 
of this model generates the covering surfaces with the corresponding combinatorial factors. 
These discrete surfaces can be interpreted as cell decompositions of continuum surfaces 
covering the target space. It should be intuitively clear that the combinatorics of these 
surfaces should not depend on the cellular decomposition of the target manifold, but only 
on global features like the topology of the manifold, the topology of the covering Riemann 
surface, and the number and types of branch points allowed. Our solution below will 
confirm this intuition. 

The solution of the model is given in terms of a sum over polynomial representations of 
the group U(N). This sum is, by construction, a generating function for the combinatorics 
of covering maps. The order of the representation gives the degree of the (connected and 
disconnected) coverings. The sum depends on the genus G of the target manifold, and 
a number of variables tied to geometrical data of the covering maps: N~ 2 is the genus 

expansion parameter of the covering surfaces, and we associate, for each cell corner p, a 

(p) 

weight r corresponding to a branch point at p of order k. 

For some applications, one would consider a slightly more general problem of counting 
coverings with branch points that can occur anywhere on the smooth target manifold. In 
order to allow for this possibility, we are led to take a continuum limit: We simply cover the 
target manifold by a microscopically small cell decomposition, with the weights of these 
branch points tuned correspondingly. We will find that the combinatorics of the resulting 
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statistics of movable branch points actually simplifies significantly over the general case of 
fixed branch points: Many special configurations of enhanced symmetry, corresponding to 
coalescing branch points, are scaled away. The weights associated with the punctures are 
not subdued to scaling; this makes the difference between the punctures and the rest of 
the points of the cell decomposition. 

Apart from the obvious mathematical interest of our approach, we are able to connect 
our results to recent work on the QCD string in two dimensions. In the case of a target 
space with nonnegative global curvature, the chiral (i.e., orientation preserving) sector of 
the Gross and Taylor || string theory describing two dimensional Yang-Mills theory, is a 
particular case of the string theory defined by our matrix model. We discuss the problem 
of the "O factors" in the string interpretation of the Yang-Mills theory and interpret the 
"O-points" as punctures in the target space. 

2. Definition of the model 

Consider a smooth, two dimensional, compact, closed and orientable manifold M. G 
of genus G with Ao marked points (punctures), which we denote by p = 1,...,Aq. The 
manifold is compact in the following sense: We assume that there is a volume form dA on 
M.G and the total area At = J dA is finite. We will consider the ensemble of nonfolding 
surfaces covering M.g and allowed to have branch points at the punctures. These surfaces 
are smooth everywhere on M.q and are given a volume form inherited from the embedding. 
In this way, the area of a surface covering n times the target manifold is equal to uAt- 

We will resolve the problem by discretizing the target manifold. Introduce a cell 
decomposition of the original target manifild M.g such that each cell is a polyhedron 
homeomorphic to a disc. The vertices of the cell decomposition are by construction the 
Ao punctures of M.g- The resulting polyhedral surface (cellular complex) Mq is thus 
characterized by its genus G, and by its set of points p, links £ and cells c. The numbers 
of points, links and cells which we denote correspondingly by Ao, Ai and A2, are related 
by the Euler formula 

Ao - M + A2 = 2 - 2G. (2.1) 

Each cell contributes a fraction A c to the total area At of the target manifold, so that 
At — Ylc A given manifold can be discretized in many different ways, but the choice 
of discretization is irrelevant for our problem. 
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A branched covering of Mq of degree n is a surface E obtained by taking n copies of 
each of the polygons of Mq and identifying pairwise the edges of the n polygons on either 
side of each link. The Riemann surface obtained in this way can have branch points of 
order k (k = 1, 2, n) representing cyclic contractions of edges. They are located at the 
points p € Mq. The discretized surface E has nH\ links, 71A/2 polygons and uMq — J2 P ^p 
points (with b p the winding number minus one at point p). Its total area is nA?- Its genus 
g is given by the Riemann- Hurwitz formula: 

2g-2 = n{2G-2) + ^b p . (2.2) 

v 

The partition function is defined as the sum over all possible coverings E — > Mg conserving 
the orientation. A factor N 2 ~ 2g is assigned to the genus g of the covering surface. Fur- 
thermore, we introduce Boltzmann weights associated with its branch points. The weight 

of a branch point of order k is t^', where k > 2. A regular (analytic) point gets a weight 

Ap) 
i 1 . 

The partition sum is now defined as a sum over all (not necessarily connected) cover- 
ings E — > Mq: 

Z= e- nAT N 2 - 29 Y[H(t^) nk . (2.3) 

S^M G p=lfc>l 

where, associated with the point p G Mq, nk{p) (with k > 2) is the number of the branch 
points of order k of E and ni(p) the number of regular points. The symmetry factor of 
the map is understood in the sum. 

Now we introduce a matrix model whose perturbative expansion coincides with ( |2.3|) . 
To each link I = (pp') we associate a field variable $^ representing an iV x N matrix 
with complex elements. By definition $ <pp / > = ^ > < p / p> - In order to be able to associate 
arbitrary weights to the branch points we will associate an external matrix field with the 
corners of the cells. Let us denote by (c,p > the corner of the cell c associated with the 
point p. The corresponding matrix will be denoted by Br c p> . 
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Fig.l: The matrices associated with the links < 12 >, < 13 >, < 14 > 
and the corners (214 >, (312 >, (413 > associated with the point 1. 

The partition function of the matrix model is defined as 

Z = I Y[[V$ e } 11 exp(e-^ N Tr$ c ), (2.4) 

J t c 

where $ c denotes the ordered product of link and corner variables along the oriented 
boundary dc of the cell c 

* c = n (2-5) 

e, P edc 

and the integration over the link variables is performed with the Gaussian measure 

JV 

[V$i] = (N/n) N2 J] d^d^ e~ N Tr (2.6) 
The perturbative expansion of fl2.4j ) gives exactly the partition function ( |2.3j ) of 



branched surfaces covering Mq. The weight tl of a branch point of order k (k > 2) 
or regular point (k = 1) associated with the vertex p e Mq equals 



tl P> = ^Tr[S p fc ], (2.7) 
where we have defined the matrices B p as the ordered product 

B v = W B {c, P > ( 2 - 8 ) 

c 

of the S-matrices around the vertex p. 
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3. Exact solution by the character expansion method 

The method consists in replacing the integration over complex matrices by a sum over 
polynomial representations of U(N). Applying to ( |2.4j ) the same strategy as in ref. @, we 
expand the exponential of the action for each cell c as a sum over the Weyl characters x h 
of these representations: 

exp(e-^ N Tr$ c ) = £ X h (<&c) e~ A ^. (3.1) 

h 

The representations are parametrized by the shifted weights h = {hi, hi-, ■ • • , Jin}, where 
hi are related to the lengths mi, of the rows of the Young tableau by hi = N — i + rrii 
and are therefore subjected to the constraint hi > h% > . . . > hjsr > 0. We will denote 
by \h\ = T^nii the total number of boxes of the Young tableau. The dimension of the 
representation h is given by 

A, = n^^ (3-2) 
- LJ - ? — % 

i<j 

and the "Omega factor" Vth by 

1=1 

An explicit representation of the Weyl characters \h * s only needed for the derivation of 
two essential integration formulas (fission and fusion rule, respectively): 

J [V*] Xh (*i $ $ 2 $*) = Xh ^ t ) X J$ 2 ), 
J [2?$] x,($i *) $2) = 6 hth . ^ x h ($i $2). 



(3.4) 



For a simple proof of these identities, see |]T|. It is now possible to exactly perform the 
integration with respect to the gaussian measure (|2.6| ) at each link. Using the fusion rule 
one progressively eliminates links between adjoining cells until one is left with a single cell. 
The remaining links along that plaquette are integrated out by applying the fission rule. A 
similar procedure was first used in the exact solution of two-dimensional Yang-Mills theory 
0,|J. Employing the very useful relation 



one finds the exact solution of the matrix model 

2-2G A/o 



h ^ ' p=l 



v) 



x h {M 



e-^ AT . (3.6) 



The characters x h (-Bp) are related to the branch point weights t£ at the site p through 
the Frobenius formula (we omit the index p for clarity): 



|/i|! _/v ni+7t2 "'~ r! ' 3 " , ~' 
ini ni i 2^ n2 ! 3n 37l3 ! 



ni+2ri 2 +3n 3 ... = |/!,| 

(3.7) 

The sum is over all partitions of the covering number \h\, and the ch/ l (l ni , 2 712 , 3 713 , . . .) 
are the characters of the symmetric group S\h\ corresponding to the representation h and 
the partition class with cycle structure (l ni , 2 n2 , 3 n3 , . . .) . Let us emphasize that (|3.6|), 
together with (|3.7|) , is the complete solution of the combinatorial problem of counting 
branched covers 



In appendix A, we have given the first few terms of the expansion (|3.6|) , in the special 
case where the branch point weights are identical on all Ao sites. We have also given the 
first few terms of the free energy F — log Z, corresponding to connected surfaces. 

A special case of ( |3.6| ) is obtained when no branch points are allowed on the target 
manifold. This is done by choosing all matrices B p = A\. Then the solution ( |3.6j ) becomes 
completely independent of the cell decomposition of the manifold: The cell corners are 
indistinguishable from any other point on the surface. We obtain 

Z = Y (^) 2G 2 e~^ A -. (3.8) 

This is completely trivial for the sphere and very simple (see 0,0) for the torus, but 
highly non-trivial for higher target space genus G > 2: 

(N 2 e~ AT , forG = 0; 

log Z = A° EZi,m=i I z~ kmAT for G = 1; (3.9) 

[ N 2 ~ 2G e~ AT + (N 2 ~ 2G ) 2 |(4 G - l)e" 2 ^ + . . . for G > 2. 

For G > 2, logZ counts the number of smooth and locally invertible maps of a surface of 
genus g onto a surface of genus G. It would be interesting to find an analytic expression 
for logZ for all G > 2. 
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4. Continuum limit: statistics of movable branch points 

In the preceding section we have presented the full solution to the problem of counting 
the branched covers of a smooth, closed and orientable target manifold of any genus G 
and Mo punctures. By construction, the solution ( |3.6|) depended, aside from G and the 
genus g of the branched cover, on the number and types of branch points at the punctures 
of M.g'- At each puncture p we are to specify a set of weig hts 4 P) 

associated with the 

winding numbers of the covering surface. It is natural to consider a related, but different 
combinatorial problem: Allow the covering surfaces to have branch points anywhere on the 
target manifold. In this case we have to sum over the positions of the branch points with 
an appropriate integral measure. 

The solution of this problem is actually contained in the solution of the previous one. 
Since branching in our model is constrained to occur at the sites of the cell decomposition, 
we are led to consider a continuum limit: The cell decomposition has to densely cover the 
manifold. The simplest choice is to assign identical weights at each branching sites, that 
is, choose B p = B everywhere, and assume that all cells have the same area A c = 4? 1 . Set 



and take A/q — * oo, while holding At and the continuum couplings fixed. Thus the 
probablity for branching at a specific site p goes to zero in a prescribed way. In this 
continuum limit the configuration space of the branch points becomes infinite and their 
statistics drastically simplifies due to the fact that the probability to have more than one 
branch point associated with the same point of Mq tends to zero. Mathematically, this 
simplification is immediately seen from the partition function ( |3.6| ), and the Frobenius 
formula ( p.7|) . The product of quotients of characters exponentiates, and (|3.6|) becomes: 



t k =—TrB k = -f for k > 2, 




(4.1) 




(4.2) 



Here, in view of ( |3.7| ), the tableau dependent numbers 



are given by 



\h\\ ch /l (ll /l l- fc ,fc 1 ) 
k(\h\ - k)\ ch^lN) 



(4.3) 
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An explicit formula is 



1 n n , v 

^-^-i)...(^-Hi)n i-o- 



(4.4) 



which is quickly derived from ( |3.5| ),( P~T| ), and the identity (coming from the Schur definition 
of the characters, see Iff) 







x-AB) 



~ k 

with hi 



hi Skik 



(4.5) 



The are actually symmetric polynomials of degree k in the weights hi (see appendix B, 
where we have also listed a few examples). They are iV independent as long as \h\ < N. 

In appendix B, we have given the first few terms of the expansion (|4.2j ), as well as the 
first few terms of the free energy T — log Z, corresponding to connected surfaces. 

Finally, let us mention that we can keep the weights associated with a given number 
of points unsealed. Then we obtain the partition function of the ensemble of branched 
coverings of a target manifold of area At, genus G and Mq punctures. The evident gener- 
alization of eqs. ( |3.6|) and ( |4.2|) is 



y (- 



2-2G No 



n 



X h (B P ) 
X h (Ai) 



exp 



\-kth 



L /c=2 



\h\An 



(4.6) 



5. Relation to chiral 2D Yang-Mills theory 

The solution of the map counting problem with movable branch points considered in 
the last section has an interesting physical interpretation. Let us restrict ourselves to the 
case of only simple branch points: = for k > 3. Then, using the explicit result for £3 
(see appendix B), and choosing 

-r 2 = - A T , (5.1) 
we can rewrite the partition function ( |4.2j ) as 

where is the second Casimir of the group U(N). This is very nearly the partition func- 
tion of two-dimensional U(N) Yang-Mills theory. This theory was shown to be solvable 
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on any target manifold in [|J,||. Recently, it was demonstrated by Gross and Taylor 0, 
that YM2 can be interpreted as a string theory. This was achieved by the tour de force 
approach of expanding the exact solution in j? and interpreting the terms as string maps. 
Here we are inverting the Gross- Taylor program: We start with a theory whose interpreta- 
tion as a string theory generating covering maps from a worldsheet to the target manifold is 
manifest, and aim to derive YM2. There remain, however, two subtle differences between 
(0) and YM 2 : 

(1) In U(N) YM<i the sum over polynomial representations h is extended to a sum 
over coupled (non-polynomial) representations. This difference is easy to understand: The 
missing representations clearly correspond to orientation reversing maps, which have been 
eliminated from the start by our chiral matrix model. Indeed it became evident from the 
work of 0, that YM2 factorizes into two rather weakly interacting chiral sectors. Our 
class of models furnishes a precise realization of one such sector. 

(2) Our partition function contains, for genus G 7^ 1, the extra factor 0^ G_2 , absent 
in YM.2- This factor eliminates the so-called "O points" and "O -1 points" introduced 
originally by Gross and Taylor in order to sustain the string picture for G 7^ 1. For G = 0, 
we can add two punctures and associate weight 1 with the windings around them (i.e. take 
a matrix source B = 1). The answer is given by ( |4.6| ) with Ao = 2 and B\ = B 2 = 1 
and coincide with the partition function of the Yang-Mills theory on the sphere. The two 
punctures are the two "O points" . In summary, in the case of a target space representing 
a sphere with two punctures our theory is exactly equivalent to chiral YM2. For G > 2 
however, it is not possible to eliminate the extra factor ( the "O -1 points") by adding 
punctures. On the other hand, the "O -1 factors" have been given an interesting re- 
interpretation in the work of Cordes, Moore and Ramgoolam 0. These authors showed 
that the YM2 partition function could be rewritten as a sum over covering maps without 
non-movable, special singularities, but weighted instead with special topological invariants 
(the "Euler character of Hurwitz moduli space"). It would be interesting if this result 
could be reproduced by a simple matrix model, similar in spirit to (|2.4| ) and possessing an 
equally clear surface interpretation. 

1 Adding a puncture to the target space of the Yang-Mills theory does not change anything. 
For example, the gauge theory defined on a cylinder with Dirichlet boundaries is identical to the 
gauge theory on the sphere. 
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6. Sphere-to-sphere maps: saddle point analysis 

In the case of a target space of spherical topology G = we can apply a saddle point 
technique to the sum over representations to extract the contribution from coverings with 
spherical topology (i.e. g = 0). This was explained in detail in and leads in the 

case of the model with fixed branch points ( |3.6|) to the general Riemann-Hilbert problem 
discussed in ||. Fortunately, the case of movable branch points (|4.2| ) is much simpler. 
Introducing a continuous coordinate h = j^hi and a density p(h) = di/dhi, one finds, with 
the help of (U), 

£ dk 'h^h/ = iog{k ~ b)+ \ At + ¥' (h) > (fu) 

where the effective potential V'(h) is, for a finite number of non-zero r n 's, a polynomial 
in h whose coefficients are selfconsistently dependent on the first few moments H n = 
iV -1-n h™ of the resolvent H(h) = dh' ■ It is found as follows. Define the 
potential 

oo 1 

V{h) = J2z T kh k . (6.2) 



k 

k=2 

Then 

nh ^^ijjj^ v{se ~" ( ' ,) - (6 ' 3) 

Here the contour surrounds the cut [6, a] of H{h). The simplest non-trivial example consists 
of taking = for k > 3. The analysis then becomes very similar to the one for chiral 
YM2, obtained in ||. Eq.( |6.1|) then reads: 

dh'p(h') , At T2,_ tX , , N 

The term log(/i — b) is a consequence^ of the fact that the hi are a set of ordered integers. 
The density is therefore constrained to have a maximum value of 1 (see H). This equation 
is solved in the standard way by a contour integral 
dh'p(h') 



H(h) 



h-h> 

h \ m ^71 — nl ds ln(s-b) + ^f + f(l-h) 

(h- s )y/(s-a)(s=b) 



= \n( T ^- b ) + v /(h-a)(h-b)f 



(6.5) 



2 A brief argument shows that for the phase corresponding to the sum over surfaces, a part of 
the density, starting at the origin and finishing at a point b, is saturated at it's maximum value. 
The sum over coverings is an expansion in powers of e~ Ar and T2. For small T2 and small e~ AT 
i.e. large At the e At ^-^ hl term in ( |4.2| ) attracts all the hi towards the origin saturating the 
constraint hi + \ < hi, p(h) < 1. 
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Expanding the contour out to infinity we catch a pole at s = h, the discontinuity across 
the cut of the logarithm, and a contribution from the contour at infinity. The final result 
is 



H(h)=lnh + ^-^ + ^(h-V(h-a)(h-b)) 



- 2 In 



\Jh — a + \Jh — b 



(6.6) 



\/a — b 

The density p(h) is as given by the discontinuity of H{h) across its cut: i 

p{h) = ^cos" 1 (J^\ - jV(h-a)(h-b), (6.7) 
and the cut points a and b are determined from the behaviour of H(h) for large h: 



Imposing this asymptotic behaviour on ( |6.6| ) leads to the two equations 

x = rle~ AT e x and r] - T 2 e~ AT+x = 1. (6.9) 



where we have defined x — At + 2 log((a — b)/A) and 77 = (a + b)/2. Differentiating 
w.r.t. to At leads to -§^F = — < h > +| where the free energy T is defined here by 
T = 1/A^ 2 ln^. The expectation value < h >= J Q a dhp(h)h can be calculated from the 
expansion for H(h) (|6.6|) . Using ( |6.9|) it is possible to integrate up to obtain an explicit 
expression for the spherical contribution to T 

^ = e-^+x(l-? x + i x a ). (6.10) 

^From (|6.9|) it is clear that x is a power series in r|e _j4T . We can thus perform a standard 
Lagrange inversion on (|6.9|) ,( |BTT0"D , and obtain after a short calculation: 



n n 3 



^=J2 ^- r ri n - 2 e~ nAT . (6.11) 



This result was first obtained in 0,0. To discuss the convergence properties of ( |6.11| ), 
it is natural to take the continuum coupling proportional to At, since the branch points 
can be located anywhere on the manifold (see also flOl) ): T2 = tAx- Note that the series 
is only convergent for t 2 A T e~ Ar < e _1 . Beyond this point the boundary conditions 



Specifically p(h) = 4-{H{h + it) — H{h — it)) where t is a small positive real number. 
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lead to a non-physical complex value for \- We see that the sum over branched coverings 
is convergent for both large and small areas. For large enough t and intermediate values 
of the area At, however, the entropy of the branch points is sufficient to cause the sum 
to diverge. It is interesting to understand this divergence in terms of the Young tableau 
density p(h). Along the critical line (T2, At) and T2 > the density becomes flat at its 
upper end point a, i.e. the singularity at the end point changes from p(h) ~ (a — h) 1 / 2 
to p(h) ~ (a — /1) 3 / 2 . Along the critical line (t2,At) and r 2 < it is the singularity of 
the density at the point b that changes from 1/2 to 3/2. This is just as occurs in matrix 
models of pure 2D gravity, and indeed for t 2 A 2 e~ AT ~ e _1 the free energy ( |6.11| ) behaves 
as 

F ~{e~ x -t 2 A T e- AT f 12 . (6.12) 

So far in our analysis we have ignored the constraint that we are summing over positive 
Young tableaux i.e. that b > 0. Since b = is not a singular point in the boundary 
conditions ( |6.9|) it does not correspond to a singularity in the sum over surfaces. If, 
however, we take the sum over representations ( |4.2| ) as fundamental (as would be the case 
for QCD2) then one should take this constraint into account. 




Fig 2. Phase diagram in the t, At plane. 
A typical density, p(h), is sketched in each phase. 
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Setting b = in ( |6.9|) leads to the pair of equations 

f 2 A 2 e - a t = xe - x and At = x + 2 i n (2 _ v /^). 



(6.13) 



which determine t and At parametrically in terms of \. This curve is plotted in Fig. 2. 
It connects the dot on the left to the shaded area on the right. Immediately below this 
line the support of the density is entirely positive, i.e. it starts at a positive value of h. In 
addition it is less than one for the full range of its support, see Fig. 2. The phase transition 
across this line is exactly analogous to the large iV phase transition that occurs in QCD 
(see [ 10l||ll|| ). There there is a phase transition separating a strong coupling regime from 
a weak coupling phase, and as is the case here, there is no singularity in either phase 
indicating the the transition point. Indeed it is possible to formulate the Gross- Witten 
||1C|| and Brezin-Gross [3 1| models in terms of a sum over representations and the large 
N phase transition is precisely the point at which the density just begins to touch (or 
pull away from) the origin. Further analysis shows that there are two further phases for 
the model, as first observed in [||. There is one phase where the density has an entirely 
positive support but attains its maximal value p(h) = 1 over a single finite interval, and 
another phase where the density starts at the origin and has two separated intervals where 
p{h) = 1. Both of these phases involve two separated nontrivial cuts and can be calculated 
in terms of elliptic functions. We do not present here explicit results for any of these 
extra phases. The complete phase diagram for the partition function is shown in Fig. 2. 
In each phase is sketched a typical density. Along the transition lines are indicated the 
corresponding critical behaviours of the density. The exact position of the line along which 
the point d equals zero has not been calculated. We indicate this by the using a dotted 
line. The sum over representations is thus seen to have a much richer phase structure than 
the perturbative sum over surfaces it generates. 



7. Concluding remarks 

We have introduced a new class of matrix gauge models which solve the general prob- 
lem of counting branched covers of orientable two-dimensional manifolds with specified 
branch point structure. The result is given as a weighted sum over polynomial repre- 
sentations of U(N). We conjecture that a similar methodology could be applied to non- 
orientable surfaces by replacing the complex matrices by real matrices and the group U(N) 
by 0(N). 
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Our approach actually treats also the case of manifolds with boundaries. Indeed, due 
to the invariance with respect to area-preserving diffeomorphisms, the problem does not 
depend on the length of the boundaries and the latter can be considered as punctures. 

As an interesting by-product of our investigation, we have derived in a constructive 
and transparent fashion important features of the original Gross- Taylor interpretation of 
YM 2 string theory. 

An interesting problem which we did not consider in this paper is the calculation of 
Wilson loops. The Wilson loop functional defined as the sum of all nonfolding branched 
surfaces spanning a closed contour C G M.q should satisfy a set of loop equations with a 



contact term similar to those considered in [13 
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Appendix A. Free energies for Aq fixed, identical branch points 



In this appendix we give some examples of how to extract the combinatorics of con- 
nected coverings of a smooth manifold with Ao fixed, identical branch points. The partition 
sum ( |3.6| ) for this case becomes 



( A h N 


2-2G 


\x h (B)] 






) 


IxMi). 





\h\A T 



(A.1) 



Define 



Z = l + J2 Z r, 



-tlAt 



(A.2) 



71=1 



Using the Frobenius formula for characters, the first few Z n 's, up to order four, are found 
to be: 



Zi = N 2-2G t No 



6 



6 



N 2 N 3 
3 



AT 



iV 2 z N 3 
3 

iV 2 " 



N 2 N 3 



(A.3) 



These expressions allow us to obtain explicit results for the map counting problem. The 
free energy counting connected surfaces is 



F = \ogZ = Y,F n e 



tlAt 



(A.4) 



n=l 



Let us present some explicit low order results: 
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G = 0: 



F 1 =N 2 tf > 

[*M>]-1 - 

Fi= £ ^ 2 - 2 ^( 2 j7 2 )(^) M, - 23 - 2 ^ +2 

3=0 * ~ + ~ 



+ 



+ iV° 



+ 2 ( ;^)^-^ 2 + i(^)^- 9 ,3 



F 4 =W 



+ C(A^" 2 ) 

I 6 / VI, 4/ 4V2,2/ 16V 3 / 16V 6 JJ 

+ ( 27 (m) +3 (m: 2 ))^- u ^+ 

a + (S))^- io *2 + 



4A/" -12 + 6 



+ 6 



+ 0(iV o ). 



G = l: 



^ \ 

9=2 



4 



7Af + 60 f W" 4 *2 + 9AT ^- 3 t 3 



(A.5) 



(A.6) 



+ C(A^" 4 ). 
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--N- 



--N- 



--N~ 



15 



[^o] +3 
3=4 



No 
2g-6 



^2^0+6-29^3-6 



-6 



220 -3M 



3 

5275 



( 



M + N -io 



+ 0(N- 10 ) 



3760W + 41280 



No' 



+ 



-12\ 



(A.7) 



Here have employed the standard notation for binomial and multinomial coefficients. 
It is instructive to draw the Riemann surfaces corresponding to the various terms in 
( |A.5|) , QO)D , (|A.7|) ■ As will be seen, the combinatorics involved increases rapidly in com- 
plexity 

The above examples are easily checked to be in agreement with the Riemann-Hurwitz 
formula. Of course, this formula only gives a necessary condition for the existence of a 
Riemann surface. The above results can be used to decide whether a Riemann surface of 
given G,g and branch point structure t^ 2 t^ 3 . . . actually exists. For example, we see from 
F4 in (|A.5| ) that there exists a fourfold cover of the sphere by a sphere with exactly six 
simple branch points, where two branchpoints are located at each of three locations of the 
target manifold (the associated symmetry factor is A). As a second example, we see from 
F% in (|A.6|) that there exists a triple cover of the torus by a double torus with exactly one 
branch point of order 2 (the associated symmetry factor is 3). 



Appendix B. Free energies for an arbitrary number of movable branch points 

In this appendix we give some concrete examples of how the maps are counted in the 
continuum limit we defined above. Define 

00 

Z = l + Y,Zn e~ nAT . (B.l) 

n=l 
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The first few auxiliary (see (|4.4|) ) are: 



iV(iV — l)(2iV — 1) 
6 



(B.2) 



3 

1 ,9 



+ T ,^-g JV+2)Eft< - jV ( JV - 1 )( 3jv2 - 3jV + 2 ) 



3 X 2 2 

They are easily computed from the following representation of eq.( 

= t? f ti h{h - l] ■ ■ ■ {h - n + x) exp li^i^r'm)] 



(B.3) 



where we have introduced H(h) = 5^£Li h \. ■ The first few Z n } s, up to order four, are: 
Z x =N 2 ~ 2G 



Z 2 =(-N 2 ) 2 ~ 2G (e^ +e -^) 



Z 3 =(-N 3 ) 2 ~ 2G ( e % T *+^ rs 
6 

+ (IiV 3 ) 2 - 2G e-** 7 * 
V 3 ; 



^2" '3 -J- e -W T 2 + 7h T 3 



)+ 



(B.4) 



1 



£ 4 = (^iv 4 ) 2 - 2G ( e ^ T2 +l^ T3 +^ T4 



'24 

+ Nr~ 2G ( 



)+ 



+ 



1 

12 



iV 4 ) 2 " 2G e -i^ T3 . 



This gives the following continuum free energies (as can be checked easily from the results 
of appendix A): 
G = 0: 

Tx =N 2 



2g+2 



6 



+ 



(B.5) 



+ 



1 



1 



+ 0(iV o ) 



1 .6 , 27 , 



-2^2 



(T 2 + 24 T2T3 + 2 T2T " 3 + T2T3T4 + 6 Ts + + 3 T2T4 + 8 T4 



+ 0(iV o ). 
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G=l: 



G = 2: 



ft 



ft 



9=2 



2g-2 



{2g-2)\ 2 
ft =| + N~ 2 [srf + 3r 3 ] + 0(iV" 4 ) 
•F 4 =^ + N~ 2 [S0ri + 9r 3 ] + 0(N~ 4 ). 



ft =AT 



15 



2-2g_ 



2g-Q 



g=4 



(2g-6)\ 



ft =N 



-6 



220 



+ iV" 8 320r 2 2 + 135r 3 + OiN' 10 ) 



J- 4=iV -8^ + iV-10 



20640r 2 2 + 8505r 3 



+ 0(N~ 12 ) 



(B.6) 



(B.7) 
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